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In the present paper we report comprehensive analysis of the thermodynamic properties of novel 
hydrogen dominant potassium hydride superconductor (KHe). Our computations are conducted 
within the Eliashberg theory which yields quantitative estimations of the most important ther¬ 
modynamic properties of superconducting phase. In particular, we observe, that together with 
the increasing pressure all the thermodynamic properties decrease, e.g. Tc € (72.91,55.50) K for 
p G (166,300) GPa. It is predicted that such decreasing behavior corresponds to the decreas¬ 
ing hydrogen lattice molecularization with increasing pressure value. Futhermore, by calculating 
the dimensionless thermodynamic ratios, familiar in the Bardeen-Cooper-Schrieffer theory (BCS), 
it is proved that KHe material is a strong-coupling superconductor and cannot be quantitatively 
described within the BCS theory. 
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I. INTRODUCTION 

In the terms of research on the phonon-mediated su¬ 
perconductivity, the hydrogen-based materials [T] expe¬ 
rience at present the apogee of the scientific attention [2] . 
The advent of this high interest is marked by the theo¬ 
retical work of Ashcroft [3] who predicted that metallic 
hydrogen, under the influence of high-pressure, can be a 
superconductor characterized by the high transition tem¬ 
perature (Tc). Such desirable character of the supercon¬ 
ducting phase in hydrogen stems from its light atomic 
nuclei and the absence of the inner electron shells, re¬ 
sulting in the high value of the electron-phonon coupling 
constant (A). However, the pressure value required to 
induce the superconducting phase in pure hydrogen is 
beyond any commercial applications. 

In this context, the introduction of the heavier atoms 
into the crystal lattice allows to lower the metallization 
pressure, preserving at the same time the extraordinary 
superconducting properties of hydrogen. In this context, 
it is possible to supplement the high-temperature super¬ 
conductors [4], by the ones described within the well- 
established electron-phonon scenario 0, 0- The true 
testimonial of the importance of the hydrogen-based su¬ 
perconductors are recent experimental results for the hy¬ 
drogen sulfide superconductor (H 2 S), which report the 
extremely high transition temperature equal to ~ 190 
K, obtained at the relatively moderate pressure value (p) 
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0: i- If reproducible, these results may constitute a 
significant breakthrough in the research on the room- 
temperature superconductivity. 

One of the newest attempts in the research on the 
superconducting properties of the hydrogen-based com¬ 
pounds consider also application of the alkali metal 
dopants. The most recent theoretical approach concen¬ 
trates on using potassium, leading to the thermodynami¬ 
cally stable hydrogen-rich layered material where each six 
hydrogens hold one potassium atom (the KHg compound 
within the C2/c phase) [S]. In particular, it is theoreti¬ 
cally predicted that superconducting phase, in this com¬ 
pound, is stable from ^ 166 GPa up to around 300 GPa. 
Please note that the upper pressure value constitutes the 
accuracy limit for the projector-augmented plane-wave 
potentials m used within the VASP code m for geome¬ 
try optimization and electronic calculations in [9] . Below 
the metallization pressure, but above 70 GPa, the KHg 
crystallizes in the metallic C2/m phase. In what follows, 
the synthesis pressure for KHg is expected to be much 
lower then in the case of the other recently predicted al¬ 
kali metal hydrides, namely the LiH-family [T2]. This 
fact arises from the faster transfer of electrons to the hy¬ 
drogen in the KHg when comparing to the LiH-hydrides, 
and it is an undisputed advantage of this compound. 

In this spirit, present paper reports the analysis of the 
thermodynamic properties of the KHg superconductor, 
which should be of crucial importance and interest for 
further design of possible hydride superconductors. The 
calculations are conducted here for three selected pres¬ 
sure values (166 GPa, 230 GPa, and 300 GPa) which 
sample the entire known superconducting phase of the 
KHg compound. Due to the relatively high values of 
the electron-phonon coupling constant for all considered 
pressure values, computations are conducted within the 
strong-coupling generalization of the Bardeen-Cooper- 
Schrieffer theory (BCS), namely in the framework of the 
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Eliashberg equations [T3]. This theoretical model allows 
us to provide the quantitative estimations of all of the 
most important thermodynamic properties of the super¬ 
conducting phase such as critical temperature, band gap 
at the Fermi level, specific heat, thermodynamic criti¬ 
cal field, and the electron effective mass. We note that 
our analysis is based on the isotropic Elishaberg equa¬ 
tions, following the character of the Eliashberg function 
{a‘^F{u})) presented in [S], and adopted for calculations 
in the present work. 


II. THEORETICAL MODEL AND 
COMPUTATIONAL DETAILS 

In order to estimate all desirable thermodynamic prop¬ 
erties of the discussed KHg superconductor we solve the 
Eliashberg equations on the imaginary axis and in the 
mixed representation by using the iterative method pre¬ 
sented in [14]. 

As already mentioned, our calculations based on the 
Eliashberg spectral functions calculated in [5], by us¬ 
ing the QUANTUM ESPRESSO code [15] with the 
norm-conserving pseudopotentials within the general¬ 
ized gradient approximation. Such spectral functions 
yields the electron-phonon coupling constants equal to: 
Ap=166GPa = 0.92, Ap=230GPa = 0.89, and Ap=300GPa = 
0.79. Futhermore, we assume that the electron deparing 
correlations are described in the terms of the Coulomb 
pseudopotential (^) equals to 0.1 for all considered cases, 
and that the cutoff frequency is set as uic = lOUmaxj 
where Umax stands for the maximum phonon frequency 
equal to: 407.76 meV at 166 GPa and 411.92 meV at 230 
GPa and 300 GPa. 

To assure the required precision, our computations are 
conducted for the 2201 Matsubara frequencies: Wm = 
^(2m — 1), with P = l/ksT, and fcs denoting the 
Boltzmann constant. In what follows, all thermody¬ 
namic properties of interest are described quantitatively 
for T > To = 5 K. 


III. NUMERICAL RESULTS 

In Fig. [^we present obtained imaginary-axis results of 
the (A) maximum value of the order parameter (Am=i), 
(B) the wave function renormalization factor (Zm=i), 
and (C) the normalized electron effective mass {m*jme, 
where m* ~ Zm=i'me is the electron effective mass, and 
rrie denotes the bare electron mass). 

In particular, for all considered pressure cases we ob¬ 
serve decrease in the values of the Am=i, Zm=l^ and 
m*/TOe functions with the increasing pressure. In what 
follows, such a decrease is also observed in the terms of 
the superconducting transition temperature calculated 
by using the relation: Am=i('7c') = 0, which gives the 
Tc values equal to 72.91 K, 70.73 K, and 55.50 K at 166 
GPa, 230 GPa, and 300 GPa, respectively. We also note, 




FIG. 1: The maximum values of the order parameter (A), 
wave function renormalization factor (B) and the normal¬ 
ized electron effective mass (C). First two parameters are 
presented as a function of temperature for selected values of 
pressure. The latter observable is plotted against the pressure 
value. 


that this decreasing behavior of Tc is slightly more ev¬ 
ident in the case of our results then in the predictions 
presented in [9], which based on the McMillan formula 
m- Furthermore, results presented in [3] show that for 
higher electron-phonon coupling constants McMillan for¬ 
mula underestimate the Tc values in KHg, whereas in the 
case when A is close to the weak-coupling limit (A < 0.3 
[HI) the overestimation can be noticed. This is due to 
the fact that BGS theory omits the strong-coupling ef¬ 
fects, which are present in the Eliashberg theory. 

The corresponding results for the value of energy gap 
at the Fermi level are presented in Fig. as a func¬ 
tion of pressure. We note that, in the first approxima¬ 
tion such calculations can be done on the basis of the 
imaginary-axis results presented in Fig. as 2Aj,^=i(0), 
where Am=i (0) ~ ^m=i (To)- However, in order to ob¬ 
tain the physical value of the energy gap (2A(0), where 
A(0) ~ A (To)) it is required to analytically continue the 
imaginary-axis results on the real axis (w). In this con¬ 
text the physical value of the energy gap can be deter- 
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mined by using A (T) = Re [A (w = A (T), T)]. 



FIG. 2: The physical value of the energy gap at the Fermi 
energy (2A(0)) as a function of pressure. The inset depicts 
pressure dependance of the 2A (0) /ksTc characteristic di¬ 
mensionless ratio, which is compared with the predictions of 
the BCS theory. 

Furthermore, in the inset of Fig. we depict deter¬ 
mined values of the characteristic dimensionless ratio of 
zero-temperature energy gap at the fermi level to the crit¬ 
ical temperature (i?A = 2A (0) /ksTc). We note that for 
all considered pressure values the Ra ratios notably ex¬ 
ceed limit set by the BCS theory = 3.53), [5], [6]. 

In particular, Ra equals to 3.98, 3.94, and 3.80 for 166 
GPa, 230 GPa, and 300 GPa, respectively. 

In the next step, our calculations concentrate on the 
determination of the normalized free energy difference be¬ 
tween the superconducting and normal state (AF/p(0), 
where /9(0) is the electron density of states at the Fermi 
level). For this purpose, we use the following relation: 

AF 

Jm ~ 

X 

where and Z:^ represent the wave function renor¬ 
malization factors for the superconducting (S) and nor¬ 
mal (N) state, respectively. Obtained numerical results 
are presented in the lower panel of Fig. (A) for all 
considered pressure values. We note, that the AF/p{0) 
function takes negative values for T G {Tq,Tc). This 
fact ensures the thermodynamic stability of the super¬ 
conducting phase in the KHg material. 

By using the determined AF/p (0) function we are next 
able to estimate the value of normalized thermodynamic 
critical field {Flc/\Jp (0)) as: 
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FIG. 3: (A) The normalized free energy (AF/p(0)) and 

thermodynamic critical field {Hc/\/p (0)) as a function 
of temperature, for selected pressure values. (B) The 
TeC^ (Tc) /Hq (0) characteristic dimensionless ratio as a 
function of pressure, along with the predictions of the BGS 
theory marked by horizontal line. 


The corresponding results are presented in the upper 
panel of the Fig. (A) as a function of temperature. 
Again calculations are conducted for the three consid¬ 
ered pressure values. For both AF/p (0) and F[c/\/p (0) 
parameters their absolute value decreases together with 
the increasing pressure. The inverse behavior is observed 
for the dimensionless ratio which corresponds to the 
zero-temperature thermodynamic critical field (Rh = 
TeC^ (Tc) /Fg (0), where H/. (0) - ijg (Tq)), as pre¬ 
sented in Fig. [^(B). Moreover, similarly like in the case 
of the Ra parameter, the values of the Rh ratio differ no¬ 
tably from the BGS predictions = 0.168), [5], [5]. 

In particular, Rh = 0.149, Rh = 0.150, and Rh = 0.156, 
for 166 GPa, 230 GPa, and 300 GPa, respectively. 

In the last step, we analyze the temperature behavior 
of the specific heat for the superconducting (C'®) and 
normal (C^) state. In what follows, the specific heat of 
the normal state is calculated as: 


(T) _ 7 
kBpiO) /3’ 
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FIG. 4: (A) The specific heat of the superconducting (G'®) 
and normal (C^) state as a function of temperature for three 
considered pressure values. (B) The AG (Tc) /C^ (Tc) char¬ 
acteristic dimensionless ratio as a function of pressure. The 
BCS predictions are marked by horizontal line. 


where 7 = (1 -|- A) denotes the Sommerfeld constant. 

In this manner the specific heat of the superconducting 
state (C^) can be obtained by using the following rela¬ 
tion: 


= C^ + AC, (4) 

where AC is the difference between the specific heat of 
the normal and the superconducting state, which can be 
calculated as: 


AC{T) l d^[AF/p{0)] 

kBp{0) P dikeTf ■ ^ 

The results for and C'® functions are presented in Fig. 
|^(A). For all considered pressure values, we observe the 
characteristic specific heat jump at the critical tempera¬ 
ture, which is marked by the vertical line. Moreover, like 
in the case of previously calculated thermodynamic prop¬ 
erties the function decreases with the increasing pres¬ 


sure. Finally, determined dimensionless ratio for the spe¬ 
cific heat (i?c = AC (Tc) /C^ {Tc)) equals Rc = 1.84 
for 166 GPa and 230 GPa, and Rc = 1.62 for 300 GPa, 
respectively, as presented in Fig. |^(B). Again obtained 
results exceed limit set by the BCS theory = 1.43), 

0 , 0 . 


IV. SUMMARY 


In this paper, by using the Eliashberg formalism, all 
the most important thermodynamic properties are deter¬ 
mined for the hydrogen dominant KHg superconductor at 
high-pressure. In particular, we obtain the quantitative 
estimations of the parameters such as critical temper¬ 
ature, order parameter, energy gap at the Fermi level, 
thermodynamic critical field, and the specific heat. All 
results are given for the three considered pressure val¬ 
ues, which sample the entire known KHg superconduct¬ 
ing phase (p S (166GPa, 300GPa)). 

On the basis of our analysis we observe decrease in the 
values of the calculated thermodynamic properties with 
increasing pressure. This fact can be followed on the ex¬ 
ample of the critical temperature value which drops from 
72.91 K to 55.50 K when pressure is raised from 166 GPa 
to 300 GPa. From the physical point of view, such de¬ 
creasing behavior corresponds to the decreasing hydrogen 
molecularization and stronger electron depairing correla¬ 
tions with increasing pressure. 

Our results are additionally summarized by the cal¬ 
culations of the characteristic dimensionless ratios, fa¬ 
miliar in the BCS theory, which correspond to the zero- 
temperature energy gap at the Fermi level (Ra), the zero- 
temperature thermodynamic critical field {Rh), and the 
specific heat for the superconducting state (Rc)- We 
note, the determined values of the Ra G (3.98,3.80), 
Rh € (0.149,0.156), and Rc G (1.84,1.62) ratios at p G 
(166,300) GPa notably differ from the predictions of the 
BCS theory for the weak-coupling limit (i?®®® = 3.53, 
i?®®® = 0.168, and = 1.43). These differences 

arise due to the strong-coupling and retardation effects 
which occur in the KHg superconductor and proves that 
this material cannot be quantitatively analyzed within 
the BCS theory which lacks the proper description of 
such effects. 
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